Geometric phase kickback in a mesoscopic qubit-oscillator system 
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We illustrate a reverse Von Neumann measurement scheme in which a geometric phase induced on a quantum 
harmonic oscillator is measured using a microscopic qubit as a probe. We show how such a phase, generated by 
a cyclic evolution in the phase space of the harmonic oscillator, can be kicked back on the qubit, which plays 
the role of a quantum interferometer. We also extend our study to finite temperature dissipative Markovian 
dynamics and discuss potential implementations in micro and nano-mechanical devices coupled to an effective 
two-level system. 
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Although it is commonly believed that quantum mechan- 
ical behaviors are preclusive characteristics of microscopic 
systems, various non-classical effects have been theoretically 
predicted and experimentally observed in massive systems. In 
opto-mechanical devices IIHTI] as well as in setups involving 
superconducting qubits coupled with nano-mechanical res- 
onators Jill^, a variety of quantum behaviors are observable 
even at large temperature and in presence of strong dissipative 
processes. In general, the problem of enforcing quantumness 
under unfavorable conditions have received great attention in 
the last years, being the topic of great interest both from a 
technological and from a fundamental point of view. In this 
line of thoughts, we consider a system in which genuine non- 
classical features not only survive to temperature and dissipa- 
tion, but they are indeed induced by these environmental in- 
fluences. While quantum effects like entanglement or negative 
values of Wigner function have been treated elsewhere ifioll . 
here we focus on the study of the generation and the detection 
of a geometric phase 1 1 ill 1 211 in a harmonic oscillator. In the 
spirit of previous works ilSf . we reconsider the Von Neumann 
measurement scheme, which models the measurement process 
as a coupling between a large measurement apparatus, used as 
a probe, and the microscopic system on which the measure- 
ment is performed, under a "reverse" prospective, using the 
microscopic system (a qubit) to measure the geometric phase 
attached to the macroscopic one (a harmonic oscillator). 

To illustrate the idea we consider a general Hamiltonian 
model and defer the presentation of a physical scenario suited 
for its realization to the second part of this work. We thus 
take a two-level system with logical states {|0), |1)} (which we 
dub as the control qubit) coupled with a harmonic oscillator 
through the interaction Hamiltonian 

H = nr]\0){0\ (8> (Pe-''^ + be'f), (1) 

where b and are the annihilation and creation operators for 



the oscillator, rj is the coupling constant and i/; is a phase that 
can be externally adjusted to change the direction of the field's 
quadrature in phase space. As long as rj and i/j are constant 
in time, the unitary operator describing the conditional time 
evolution of the system is 



t/^(f) = |lKl|®l + |OKO|®D[f], 



(2) 



where DK] - e^* ^ '\ with 4" — -i C ne is a displace- 
ment operator acting on the oscillator. Eq. (|2|) describes a spin 
state-dependent displacement of the oscillator, which remains 
unperturbed if the control qubit is in |1) and is displaced by 
^ if the control qubit is in |0). By changing in time (p we 
can drag the state of the harmonic oscillator along a nontriv- 
ial path in phase space. In particular, by varying (p along a 
closed loop we can associate a purely geometric phase to the 
state of the harmonic oscillator. Following Ref. Il4 11511 . we 
assume that the value of (p is changed in n time steps 6t such 
that f = ndt, 5(i = ^,(5r, and U(t) = \\"=i D(5(i). Recalling that 
D{a)b(J3) = exp{i3(ayS*))D(Q' + p), we have 



U(t) = exp 



k=2 



(3) 



In the continuous limit we take Y!l=\ ^Ci'^C^ 
Y!k=2^^kY!iZ\ ^(.^^ J CdC and assume that tp is changed 
along a closed loop in a time r such that ^(r)=^(0). Making 
use of Stokes' theorem 111] we find U(T)^e'^D(0) with 



the area enclosed by the cyclic path in parameter space iV, 

Let us now suppose that the system can be initialized in the 
state |iA(0)> = \+}\a}, where |-H> = (1/ V2)(|0) H- 11» is the state 
of the control qubit and \a) is a coherent state of the oscillator. 
We then change (p{t) along a closed path spanned in a time r. 
The joint state of the qubit and oscillator at time r will be 



|^(r)> = (l/V2)(|l)- 



|0» ® \a}. 



(4) 
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FIG. 1: The oscillator's conditional dynamics pictured in phase 
space. In (a) the oscillator is displaced along a square whose area 
is proportional to the phase 9. In (b) the oscillator is displaced while 
undergoing a dissipative process. Here "Z/J^ and are the super- 
operators describing the unitary and dissipative evolution of duration 
6t, respectively. 



To better understand such a result, we discuss an illustrative 
example where we consider a rectangular path implemented 
by a stepwise change of (p: its value is set to zero for an inter- 
val of length T and then changed to (p„ = «7r/2 at time t - nT 
with « = 1,2,3. As shown above, the state of the qubit ac- 
quires a phase § - J[. - rf-T^ equal to the area enclosed by 
the path [see FiglT](a)]. We now stress three important points. 
First: although in the above example the phase explicitly de- 
pends on r, it is invariant with respect to the parametrization 
of the path. Indeed, if the area enclosed by the path remains 
unchanged, even for an arbitrary dependence of f(f) on time 
we would obtain the same result -& - Second: for a closed 
loop, J? does not depend on the amplitude a of the initial state 
of the harmonic oscillator Third: the phase acquired by the 
harmonic oscillator is kicked back on the qubit state, which 
takes the role of a microscopic interferometer. As a final re- 
mark, it is important to stress that the geometric phase ad- 
dressed here is a specific result of our interaction model and 
of the overall scheme that we have set up. 

We now address the first two points that have been men- 
tioned above. As shown by Eq. when the evolution of the 
system is purely unitary, the geometric phase picked up by the 
control qubit does not depend on the amplitude a of the os- 
cillator's initial state. This leaves the value of 6 unchanged 
even if the initial state of the oscillator is a thermal state 
jdaP{a,0,V)\a}{a\,wheieP{a,ao,V) - ^ -■-^'^ 



Pv 



is the Gaussian thermal distribution centered at point ao in 
the phase space and having variance V = (e''H-l)/(e''-l), 
with fi=Ha)„,/ki,T {kb is the Boltzmann constant and T is 
the temperature of the oscillator). By taking the initial state 
Po - ® Pv, it is straightforward to see that a phase 

identical to the pure-state case is acquired by the qubit. In 
fact, the state at time r reads p(t) - \ip}{(f \ ® pv, where 
\ip} = (|1) + e-''2'r^'|0))/ V2. The invariance of in such a 
mixed-state scenario has also been confirmed using the frame- 
work for the evaluation of mixed-state geometric phases pro- 
posed in Ref. llST . This shows that ignorance on the initial 
preparation of the state of the harmonic oscillator does not af- 
fect the possibility to generate and detect a geometric phase. 



We now assess the potential effects that non-unitary dy- 
namics may have on the occurrence of the geometric phase 
under scrutiny. On one hand, the consideration of an explic- 
itly open dynamics will make our proposal closer to the re- 
ality of the potential experimental situations that will be ad- 
dressed later on in this work. On the other hand, it is rea- 
sonable to expect significant deviations from the results found 
so far when we are far from unitarity. We thus consider the 
oscillator as affected by dissipation at rate y with the con- 
trol qubit still evolving unitarily, a situation that is formally 
described by the master equation p - -{i/fi)[H,p] + £.p, 
where X-p-jibph^ -{b^ b, p} /2) formally describes the oscilla- 
tor's damping. Although we do not refer to any explicit exper- 
imental configuration, at this stage, the analysis that we per- 
form here adheres very well to the experimental observations 
on the open-system dynamics of mechanical systems availbale 
to date Hi [80. In order to grasp the effects of the dissipative dy- 
namics, we divide the time-window of the evolution in small 
intervals, each of length St. Inspired by the Suzuki-Trotter for- 
mula 11191], the dynamics can then be approximated by alternat- 



ing a unitary evolution described by 'U'^^p-U^(6t)pU^(6t) and 



the purely dissipative propagator I)s,p-e ^"p SEll. After 
iterations, we have the evolved state 



p{N6tMD,r€{lfpo. 

This approach is particularly useful in analyzing a damped 
harmonic oscillator Indeed, the action of the dissipative su- 
peroperator £), on the off-diagonal elements of a density ma- 
trix written in a coherent state basis is given by the dyadic 
expression l22ll 



(5) 



(6) 



where \Aj} (j =1,2) are two coherent states. As we are inter- 
ested in short-time intervals 6t, we take 1 - exp{-y6t) ^ -ydt. 
Therefore, the action of I)^, on the state of our system re- 
sults in the displacement of the harmonic oscillator and the 
exponential decrease of its initial amplitude a. Moreover, 
from Eq. ^ we see that a phase factor is attached to the 
off-diagonal elements of the density matrix. Such features 
are useful to close the path across which the oscillator is dis- 
placed, as we now show using an example close, in spirit, to 
the previous one. Differently from the case of a unitary evo- 
lution, when dissipation is included the phase ip is set to for 
f € [toJo + Ti] and to tt for f € [to + Ti,to + Ti+ T2] ■ The state 
at time to + Ti is given by p{to + Ti) = (t)si'll%)^' Po, with 
N\6t - T\. Taking the limit 6t 0, this turns out to be 



p(fo + 7^i)=^(|l,»iXl,ail + |0,ai 



lO,^! - (^iXl,^!! +/1.C.), 



(7) 



where 



:^(1 
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r 

n = -^[yTi + -(1 

2y^ 2 



), ffi 



(8) 
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FIG. 2: (a) Probabilities and V-^ against the displacement cfo and 
the parameter V for rjly = 0.05 and yT^ = 20. (b) Same probabilities 
against V for and 0-0 = and the same parameters as in panel (a). 



Eq. (|7]i shows that the state of the oscillator at time T\ is con- 
ditionally displaced by a quantity -ifi\ . Moreover, the state 
of the qubit-oscillator system, which is damped at a rate Fi , 
aquires a phase factor 0{a) that, differently from the unitary 
case, depends linearly on a. We can then proceed to evauate 
the state of the system at time fo + Ti + by setting ip-n and 
taking p(fo + Ti + T2) = (Dsi^Ul)"'- p{to + Ti), with Nidt = T,. 
This displaces the state of the oscillator by in the oppo- 
site direction to what occurred at Ti . The time interval T2 is 
chosen such that the oscillator displacement —iPi accumulated 
during T\ is cancelled 12311 . By calling f - Ti + T2, the final 
state reads 



p(to -H f) = i [1 + e"^(e"''''"'|0)(l| + h.c.)] ® |ae"''^)<ae"''^ 

(9) 



where 0(a) -rja 



easily seen that the phase d{a) gained by the oscillator in 
this process is equal to the area ^ enclosed by the displace- 
ment path in parameter space, as shown in Fig. [T] (b). The 
detectability of such phase depends on the function F, which 
determine the decoherence rate for the off-diagonal terms 
in p(fo + r). Indeed, in order to achieve a non-vanishing 
phase, we should ensure that F<s; 1 . Remarkably, while 6{a) 
depends on the amplitude of the initial coherent state, F does 
not. Surprisingly, by choosing rjly <^ 1 (which embodies 
the weak-coupling condition between the oscillator and the 
control qubit), the requirement F <sc 1 is fulfilled. On the 
other hand, we can achieve any value of 0{a) making an 
appropriated choice for the value of a, so as to compensate 
the conditions required for a negligible damping. 

The approach described above can be applied so as to evalu- 
ate the effects, on the geometric phase, due to a thermal prepa- 
ration of the state of the oscillator undergoing dissipative dy- 
namics. We thus assume that the initial state of the oscillator is 
the displaced thermal state p^" = ^ daP{a,aa,V)\a){a\. Fol- 
lowing the lines sketched so far, we arrive at the evolved state 



Pv°(fo + r) 



^ d^aP(a, 



ao, V)p(to + f). 



(10) 



where p(fo + f) is given by Eq. (|9]). In light of the depen- 
dence of the phase 6 on the amplitude a [as shown right af- 



,(2-e-y,) ' .=/7^F(r,ri,r2)/2r2 with e 

r(j, Ti , T2) a dimensionless function that behaves as for 
y — » 0, thus ensuring that F — > as y — > l24ll . It is 



ter Eq. ([9])], the control qubit and the oscillator end up in a 
correlated state. This complicates the calculation of the over- 
all geometric phase associated with Py°(fo + t). Nevertheless, 
it is still possible to evaluate the geometric phase by adopt- 
ing the framework developed in Ref. llSll . which is based on 
the probability that a measurement over the state of the con- 
trol qubit has outcomes {|+>, |-)) (with |-> = (|0) - |1))/ V2). 
In order to understand this, let us first consider the state in 
Eq. (|9]l and suppose to project the control qubit onto the {|+)) 
basis. The corresponding outcome probabilities are given by 
V± - (±|Tr,„{p(f() + f))|+), where Tr„, denotes the partial trace 
over the oscillator's degrees of freedom. A straightforward 
calculation shows that V-^ - [1 + v cos 0{a)]l2, where v - e^^ . 
This reminds us of the fringes of an interferometer whose vis- 
ibility is v: the state of the composite qubit-oscillator sys- 
tem evolves along two branches (one associated to 1 1) and the 
other to |0)) that can be seen as two arms of a Mach-Zehnder 
interferometer. The two components of the state of the system 
that have undergone the evolution ruled by H (and possibly 
the dissipative dynamics considered here) are then let inter- 
fere by projecting the qubit state onto the superposed basis. 
This analysis offers us an operative interpretation of d(a). In 
fact, by changing it, the probabilities !P-t change, reaching the 
complete inversion when 9{a) - n. Therefore, 9{a) can be 
operatively seen as the inversion of the outcome probabilities 
We can attach an analogous meaning to the phase asso- 
ciated with a mixed initial state. We thus now consider the 
state in Eq. (fTOb . which gives us !P± = (1 + |A| cos[0])/2 with 
arg(A) and 



A 



d aP{a, ao, V)e 



ie(a) 



(11) 



As in the case of a pure state, the phase 9 is operatively 
defined through the inversion of the probabilities P^. Fig. |2] 
(a) shows such quantities against the initial displacement oq 
and the temperature V. Although the visibility of the fringes 
decreases with the increasing temperature, it is possible to see 
a population inversion even for high values of V. Another very 
interesting situation is the one in which the initial state is a 
non-displaced thermal state, i.e. the Gaussian distribution in 
Eq. ( [Tol l is centered in ^ = 0. The behavior of the outcome 
probabilities P-t against the thermal variance V is shown in 
Fig. |2] (b). A larger temperature results in an increase (de- 
crease) of ('P+)- The partial inversion of the probabilities 
is due to the average geometric phase 9(a) picked up by the 
oscillator during the process. 

In summary, we have shown how to generate a geometric 
phase on a system in which a qubit is coupled to a harmonic 
oscillator. The phase can be detected using the qubit as an 
interferometer. We propose systems combining effective two- 
level devices to mechanical modes aspotentiai scenarios for 
the implementation of our proposal |lll|9(]. In the nano-scale 
domain, Hamiltonian models of a form close to the one pro- 
posed here can be achieved by capacitively combining a nano- 
cantilever to a Cooper-pair box or growing a quantum dot on 



4 



a nano-beam ||^ 25]. At the microscopic scale, on the other 
hand, the coupling in Eq. ([T]i can be engineered by means 
of a three-level atom trapped within the volume of a pumped 
optomechanical cavity field and off -resonantly coupled to the 
latter [Iffl. Analogous configurations have been recently pro- 
posed 12611 as valid alternative to consolidated schemes for the 
coupling between a mechanical mode and the vibrational de- 
grees of freedom of a single atom, an ensemble of them or a 
levitating nanoparticle 126112711 . Under the presence of a dis- 



sipative environment and for a mixed thermal state of the os- 
cillator, the geometric phase can still be observed under con- 
ditions over the coupling between the qubit and the oscillator 
that can be matched experimentally llllsl]. For a nano-beam 
with fundamental frequency ~ lOOMHz coupled to a super- 
conducting qubit at a rate rj ~ IMHz and having a (realistic) 
decay rate of ~ 1 - lOMHz, which are values well within 
the validity of our approach, a temperature of 0.5K keeps the 
probabilities !P± at the visible level. By driving the mechani- 
cal mode with a two-tone signal 12811 . which is possible opti- 
cally and electrically, thus covering both the micro- and nano- 
scale configurations. In the first scenario, one could consider, 
for instance, a single Cs atom coupled to a light mechani- 
cal resonator (masses are typically in the ng range) in both 
the end-mirror or membrane-in-the-middle arrangement HI], 
as recently considered for the problem of coupling the ex- 
ternal degrees of freedom of an atom to the vibrations of a 
massive mechanical oscillator ll26ll . High-finesse cavities with 
small waists are currently employed in controllable optome- 
chanical experiments (finesse~10^ with a waist of a few fim), 
thus guaranteeing a strong enough light-atom interaction that 
is suitable for the achievement of the effective Hamiltonian 
model proposed in Eq. O]). For the examples discussed here, 
all the experimental observations that are currently available 
are in full agreement with a Markovian description of the dy- 
namics induced by the thermal background of phonons affect- 
ing the mechanical oscillator, thus making the our study per- 
fectly appropriate. Our proposal moves along the lines of an 
investigation assessing quantum effects in macroscopic sys- 
tems. It enlarges the fan of indicators of quantumness at the 
meso-scale with a figure of merit, the geometric phase, that 
arises in virtue of the sole coherent qubit-oscillator interaction 
and survives against plainly adverse operating conditions. 
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